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LONG-TERM  GOALS 

It  is  well  known  that  a  substantial  portion  of  the  air/sea  fluxes  of  heat,  moisture,  and  momentum  is 
accomplished  via  intermittent  processes  (Khalsa  and  Greenhut  1985),  processes  that  are  poorly 
understood  at  the  present  time.  Recently,  Mahrt  (1989)  and  Sikora  and  Young  (1993)  have 
demonstrated  that  coherent  structures  in  the  marine  boundary  layer  (MBL)  are  responsible  for  this  flux 
intermittency.  These  coherent  structure  types  include  such  secondary  circulations  as  two-dimensional 
rolls  (cloud  streets),  three-dimensional  convective  cells  (thermals),  and  shear-driven  eddies  (billows) 
(Brown  1980).  These  features  occur  in  different  atmospheric  boundary  layer  thermal  stratification  and 
shear  regimes;  some  are  forced  primarily  by  thermodynamic,  and  others  by  dynamic,  mechanisms. 

Our  ultimate  goal  is  to  determine  the  mechanisms  underlying  the  intermittency  in  air/sea  fluxes 
produced  by  these  coherent  structure  types.  As  summarized  below,  we  are  using  a  variety  of 
complementary  statistical/mathematical  approaches  to  objectively  identify  the  spatial  and  temporal 
characteristics  of  these  structures.  Our  primary  data  sources  include  both  the  high  resolution  output 
produced  by  the  Penn  State  version  of  Moeng's  Large-Eddy  Simulation  (LES)  code  (e.g.  Schumann 
and  Moeng  1991)  and  observations  from  the  MBL  ARI  experiments  performed  in  1995  off  the 
California  coast. 

OBJECTIVES 

Our  final  scientific  objective  of  the  project  was  to  complete  analysis  of  the  1995  MBL  ARI  field 
program  datasets  using  our  own  specialized  implementations  of  these  statistical/mathematical 
approaches  (Rinker  1995;  Rinker  et  al.  1995;  Rinker  and  Young  1996;  Winstead  el  al.  1995;  Rohrbach 
1996;  Suciu  1996;  Rogers  1997;  Rogers  et  al.  1997;  Shirer  et  al.  1997,  1998;  Rishel  1998;  Mason 
1999).  To  identify  the  spatial  and  temporal  behavior  of  the  coherent  structures,  we  use  obliquely 
rotated  principal  component  analysis  (PCA).  To  capture  the  contribution  of  each  coherent  structure 
type  to  intermittency,  we  use  both  the  capacity  dimension  (Takens  1981;  Henderson  and  Wells  1988) 
and  the  multiscale  line-length  algorithm  of  Higuchi  (1988);  we  provide  a  thorough  review  of  how  we 
use  these  chaos  quantities  in  the  MBL  book  chapter  by  Shirer  et  al.  (1998). 
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APPROACH 

PCA  has  been  shown  to  be  capable  of  distinguishing  and  quantitatively  describing  multivariate 
structures  within  the  atmosphere  (e.g.,  Richman  1986;  Preisendorfer  1988;  Alexander  el  al.  1993).  To 
identify  the  individual  coherent  structure  types,  we  use  both  the  spectral  coherence  of  pairs  of  temporal 
PC  score  series  as  well  as  their  chaotic  behavior.  Nonperiodic,  temporal  variation  is  chaotic  if  the 
details  of  a  particular  time  series  can  not  be  simulated  beyond  a  few  cycles  with  virtually  identical 
initial  conditions,  a  situation  typical  of  most  atmospheric  flows.  Our  hypothesis  is  that  each  coherent 
structure  type  has  its  own  identifiable  chaotic  behavior  as  revealed  by  how  the  fractal  dimension 
estimates  vary  as  a  function  of  scale.  As  we  summarize  below,  we  have  tested  this  hypothesis 
extensively  using  several  quantitative  measures  of  chaos  behavior  and  datasets  drawn  from  LES  and 
marine  surface  layer  observations. 

WORK  COMPLETED 

Our  analysis  required  that  we  both  test  existing  approaches  and  develop  new  methods  to  assess  the 
chaotic  properties  of  the  coherent  structures  provided  by  the  PCA.  We  tested  the  Higuchi  multiscale 
algorithm  (Higuchi,  1988;  Winstead  1995),  and  we  developed  two  new  algorithms  to  calculate 
estimates  of  the  fractal  dimension.  The  first  of  the  new  algorithms  estimates  the  correlation  dimension 
and  is  discussed  in  a  paper  by  Shirer  et  al.  (1997).  The  second  estimates  the  capacity  dimension;  the 
algorithm  is  developed  in  an  MS  thesis  by  Suciu  (1996)  and  summarized  in  a  book  chapter  by  Shirer  et 
al.  (1998).  Our  procedure  for  identifying  coherent  structures  was  first  tested  using  LES  data  by  MS 
students  Donald  Rinker,  Joseph  Rohrbach,  Nathaniel  Winstead  and  Laurentia  Suciu.  Our  approach  was 
next  extended  by  MS  students  Aric  Rogers,  Jeremy  Rishel  and  Richard  Mason  to  analysis  of  the 
measurements  taken  in  1995  MBL  off  the  coast  of  California.  These  coherent  structure  results  have 
been  published  in  numerous  articles  and  conference  papers  as  summarized  in  the  next  section. 

RESULTS 

Using  both  standard  and  newly  developed  PCA  algorithms,  we  studied  several  LES  datasets  to  see 
which  principal  components  (PCs),  or  coherent  structure  types,  are  largely  independent  of  the  large- 
scale  forcing  and  which  vary  sensitively  with  it  (Rinker  1995;  Rinker  et  al.  1995;  Rinker  and  Young 
1995;  Winstead  1995;  Winstead  et  al.  1995;  Rohrbach,  1996;  Suciu  1996;  Shirer  et  al.  1998).  These 
results  identified,  for  each  forcing  regime,  the  primary  physical  processes  and  coherent  structure  types 
associated  with  air-sea  flux  intermittency.  Significantly,  our  applications  of  the  PCA  algorithms  to 
idealized  data  tests  demonstrated  that  the  method  is  able  to  distinguish  multiple  coherent  structure 
types  under  several  realistic  conditions.  These  tests  provided  proof  that  PCA  can  yield  valid  results 
without  having  an  a  priori  conceptual  model,  as  required  of  previous  (conditional  sampling)  methods. 

In  the  second  half  of  the  project,  we  acquired  the  1995  RP  FLIP  MBL  ARI  dataset  and  calibration 
codes  from  Dr.  Jim  Edson  of  WHOI  and  Drs.  Carl  Friehe  and  Tihomir  Hristov  of  UC-Irvine.  Graduate 
students  Aric  Rogers  and  Jeremy  Rishel  completed  calibration  of  these  data,  which  included 
incorporating  the  new  cup  anemometer  calibration  algorithm  of  Hristov  and  Friehe  (1998)  into  the 
WHOI  calibration  code.  Rogers  (1997)  and  Rogers  et  al.  (1997)  described  a  case  study  of  the 
convective  atmospheric  surface  layer  using  these  data  and  the  coherent  structures  identification 
algorithm  described  above.  Rishel  (1998)  extended  the  work  of  Rogers  (1997)  and  showed  that 
Rogers’  single  case  study  results  are  robust.  Mason  (1999)  will  complete  the  series  of  RP  FLIP  case 
studies  by  analyzing  the  flux  contributions  of  each  coherent  structure  type  observed  in  the  experiment. 
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These  results  are  summarized  in  a  chapter  (Shirer  et  al.  1998)  of  a  book  edited  by  Gary  Geernaert  that 
summarizes  the  entire  MBL  project. 

IMPACT/APPLICATIONS 

The  observational  and  modeling  studies  performed  in  this  study  will  lead  to  improved  understanding  of 
the  flux  intermittency  commonly  observed  in  the  MBL.  This  will  help  advance  our  overall 
understanding  of  processes  that  affect  the  state  of  the  sea  surface. 

TRANSITIONS 

Improved  understanding  of  air/sea  flux  intermittency  will  lead  to  improved  means  for  interpreting  sea 
surface  roughness  patterns  on  SAR  imagery.  This  improvement  will  become  possible  as  we  obtain  an 
increased  understanding,  via  this  ARI  project,  of  the  vertical  momentum  transports  by  the  various 
coherent  structure  types.  Collaborative  work  with  Don  Thompson  and  Robert  Beal  of  JHUAPL  that 
was  funded  by  ONR  at  the  same  time  of  this  ARI  project  took  advantage  of  this  approach  to  derive 
from  SAR  imagery  quantitative  estimates  of  boundary  layer  depth,  surface  wind  speed  and  direction, 
and  air- sea  fluxes. 

RELATED  PROJECTS 

Drs.  Jim  Edson  and  Carl  Friehe  have  suggested  that  another  deployment  of  RP  FLIP  is  necessary  to 
sample  the  near-surface  regions  of  the  atmospheric  and  oceanic  boundary  layers  much  better  than  was 
possible  during  the  MBL  ARI.  The  goals  of  such  an  experiment  would  be  to  understand  the  processes 
in  the  wavy  boundary  layer  much  better  through  the  study  of  the  interactions  among  ocean  waves, 
surface  layer  oceanic  and  atmospheric  eddies,  and  mixed-layer  atmospheric  eddies.  In  our  analysis  of 
the  1995  RP  FLIP  data,  we  identified  a  coherent  structure  type  associated  with  the  interaction  between 
the  near-surface  wind  flow  and  the  ocean  waves,  but  there  were  not  enough  observations  close  to  the 
surface  for  us  to  resolve  the  structure  very  well  (Shirer  et  al.  1998).  With  observations  of  the  type 
suggested  by  Edson  and  Friehe,  our  objective  algorithm  would  be  able  to  provide  one  means  for 
quantifying  the  flux  intermittency  of  the  coherent  structure  types  in  the  wavy  boundary  layer. 

REFERENCES 

Alexander,  G.D.,  G.S.  Young,  and  D.V.  Ledvina,  1993:  Principal  component  analysis  of  vertical 
profiles  of  Q,  and  Q2  in  the  tropics.,  Mon.  Wea.  Rev.,  121,  1-13. 

Brown,  R.A.,  1980:  Longitudinal  instabilities  and  secondary  flows  in  the  planetary  boundary  layer.  A 
review.  Rev.  Geophys.  Space  Phys.,  18,  683-697. 

Henderson,  H  W.  and  R.  Wells,  1988:  Obtaining  attractor  dimensions  from  meteorological  time  series. 
Advances  in  Geophysics,  30,  205-237. 

Higuchi.  T.,  1988:  Approach  to  an  irregular  time  series  on  the  basis  of  the  fractal  theory.  Physica,  31D, 
277-283. 

Hristov,  T.  and  C.  Friehe,  1998:  Linear  time-invariant  compensation  of  cup  anemometer  inertia.  J. 
Atmos.  Oceanic  Tech ,  submitted. 


4 


Kaimal,  J.C.  and  J.A.  Businger,  1970:  Case  studies  of  a  convective  plume  and  a  dust  devil.  J.  Appl. 
Meteo  ,  9,  612-620. 

Khalsa,  S.J.S  and  G.K.  Greenhut,  1985:  Conditional  sampling  of  updrafts  and  downdrafts  in  the  marine 
atmospheric  boundary  layer.  J.  Atmos.  Sci. .  42,  2550-2562. 

Mahrt,  L.,  1989:  Intermittency  of  atmospheric  turbulence. ./.  Atmos.  Sci.,  46,  79-95. 

Mason,  R.,  1999:  Fluxes  by  coherent  structures  in  the  marine  atmospheric  surface  layer.  MS  Thesis  in 
progress,  Penn  State  University. 

Preisendorfer,  R.W.,  1988:  Principal  Component  Analysis  in  Meteorology  and  Oceanography. 
Developments  in  Atmospheric  Science,  17,  Elsevier  Press,  425  pp. 

Richman,  M.B.,  1986:  Rotation  of  principal  components.  J.  Climatol.,  6,  293-335. 

Rinker,  D.K.,  Jr.,  1995:  Use  of  obliquely  rotated  principal  component  analysis  to  identify  coherent 
structures.  MS  thesis,  Penn  State  University,  42  pp. 

Rinker,  D.K.,  Jr.,  T.D.  Sikora,  and  G.S.  Young,  1995.  Use  of  obliquely  rotated  principal  component 
analysis  to  identify  coherent  structures.  Preprints,  11th  Symposium  on  Boundary  Layers  and 
Turbulence,  Charlotte,  NC,  American  Meteorological  Society,  417-420. 

Rinker,  D.K.,  Jr.,  and  G.  S.  Young,  1996:  Use  of  obliquely  rotated  principal  component  analysis  to 
identify  coherent  structures.  Bound.  Layer  Meteo,  80,  19-47. 

Rishel,  J.,  1998:  A  methodology  for  objectively  identifying  coherent  structures  within  the  marine 
atmospheric  surface  layer.  MS  Thesis,  Penn  State  University,  51  pp. 

Rogers,  A.,  1997:  Chaotic  marine  atmospheric  boundary  layer  structures  isolated  and  identified  using 
statistical  and  temporal  analysis  techniques.  MS  Thesis,  Penn  State  University,  45  pp. 

Rogers,  A.,  H.N.  Shirer,  G.S.  Young,  L.  Suciu,  R.  Wells,  J.B.  Edson,  S.W.  Wetzel,  C.  Friehe,  T. 
Hristov,  and  S.  Miller,  1997:  Using  chaotic  behavior  of  the  time  series  observed  on  FLIP  to  identify 
MABL  coherent  structures.  Preprints,  12th  Symposium  on  Boundary  Layers  and  Turbulence, 
Vancouver,  BC,  American  Meteorological  Society,  243-244. 

Rohrbach,  1996:  The  dynamics  and  three-dimensional  structure  of  the  coherent  eddies  of  the  boundary 
layer  investigated  through  principal  component  analysis.  MS  Thesis,  Penn  State  University,  86  pp. 

Schumann,  U.  and  C.-H.  Moeng,  1991:  Plume  fluxes  in  clear  and  cloudy  convective  boundary  layers., 
J.  Atmos.  Sci.,  48,  1746-1757 

Shirer,  H.N.,  C.J.  Fosmire,  and  R.  Wells,  1997:  Estimating  the  correlation  dimension  of  atmospheric 
time  series.  J.  Atmos.  Sci.,  54,  211-229. 

Shirer,  H.N.,  G.S.  Young,  R.  Wells,  A.N.  Rogers,  J.P.  Rishel,  R.A.  Mason,  L.  Suciu,  N.S.  Winstead, 
H.W.  Henderson,  D.K.  Rinker,  Jr.,  J.W.  Rohrbach,  J.  Edson,  C.  Friehe,  S.  Wetzel,  S.  Miller,  and  T. 
Hristov'  1998:  Identifying  coherent  structures  in  the  marine  atmospheric  boundary  layer.  In  Air-Sea 


5 


Exchange — Physics,  Chemistry,  Dynamics,  and  Statistics ,  Gerald  L.  Geernaert,  Editor,  Kluwer 
Academic  Publishers,  (in  press). 

Sikora,  T.D.  and  G.S.  Young,  1993:  Observations  of  planview  flux  patterns  within  convective 
structures  of  the  marine  atmospheric  surface  layer,  Boundary  Layer  Meteorology,  65,  273-288. 

Suciu,  L.,  1996:  Estimating  the  capacity  dimension  of  time  series  produced  by  large  eddy  simulation. 
MS  Thesis,  Penn  State  University,  70  pp. 

Takens,  F.,  1981:  On  the  numerical  determination  of  the  dimension  of  an  attractor.  Lect.  Notes  Math., 
1125,  99-106. 

Winstead,  N.S.,  1995:  Diagnosing  chaotic  behavior  in  time  series  produced  by  large  eddy  simulation. 
MS  Thesis,  The  Pennsylvania  State  University,  61  pp. 

Winstead,  N.S.,  H.N.  Shirer,  H.W.  Henderson  and  R.  Wells,  1995:  Diagnosing  chaotic  behavior  in 
time  series  produced  by  large  eddy  simulation.  Preprints,  11th  Symposium  on  Boundary  Layers  and 
Turbulence,  Charlotte,  NC,  American  Meteorological  Society,  383-386. 

PUBLICATIONS 

Rogers,  A.,  H.N.  Shirer,  G.S.  Young,  L.  Suciu,  R.  Wells,  J.B.  Edson,  S.W.  Wetzel,  C.  Friehe,  T. 
Hristov,  and  S.  Miller,  1997:  Using  chaotic  behavior  of  the  time  series  observed  on  FLIP  to  identify 
MABL  coherent  structures.  Preprints,  12th  Symposium  on  Boundary  Layers  and  Turbulence, 
Vancouver,  BC,  American  Meteorological  Society,  243-244. 

Shirer,  H.N.,  C.J.  Fosmire,  and  R.  Wells,  1997:  Estimating  the  correlation  dimension  of  atmospheric 
time  series.  J.  Atmos.  Sci.,  54,  211-229. 

Shirer,  H.N.,  G.S.  Young,  R.  Wells,  A.N.  Rogers,  J.P.  Rishel,  R.A.  Mason,  L.  Suciu,  N.S.  Winstead, 
H.W.  Henderson,  D.K.  Rinker,  Jr.,  J.W.  Rohrbach,  J.  Edson,  C.  Friehe,  S.  Wetzel,  S.  Miller,  and  T. 
Hristov,  1998:  Identifying  coherent  structures  in  the  marine  atmospheric  boundary  layer.  In  Air-Sea 
Exchange — Physics,  Chemistry,  Dynamics,  and  Statistics,  Gerald  L.  Geernaert,  Editor,  Kluwer 
Academic  Publishers,  (in  press). 


